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Abstract 

The baryogenesis is reanalyzed based on the model by A.G.Cohen et 
al., in which the lepton number, generated by the neutrinos' scattering 
from the bubble walls appearing in the development of the electroweak 
phase transition, is converted to the baryon number excess through the 
sphaleron transition. A formula obtained in this paper on the lepton 
number production rate is correct for the both thin and thick walls 
within the linear approximation. 

Investigation on the time-development of the first order phase tran- 
sition is simulated, including the temporal change of the wall velocity 
as well as the fusion effect of the bubbles. The details of such phase 
transition dynamics are found to affect considerably the final value of 
the baryon number excess. 



Baryogenesis in the universe was studied originally in the grand unified 
theories in the late 70's 0. After the discovery of the baryon number and 
lepton number violation mechanism through the sphaleron transition [Q], 
we have another possibility of generating baryons at the electroweak (EW) 
scale 0. 

Among the various model of baryogenesis at the EW scale, we take up the 
model by A.G.Cohen et al. 0, in which the lepton number, generated by the 
neutrino's scattering from the bubble walls appearing in the development of 
the EW phase transition, is converted to the baryon number excess through 
the sphaleron transition. 

The purpose of his paper is to find a formula of the lepton number pro- 
duction rate which is correct for both thin and thick walls when the vacuum 
expectation value of the Higgs scaler changes linearly whitin the wall (which 
may be called as the linear approximation). 

The other aim of ours is to elucidate the time developement of the phase 
transition in detail including the temporal change of the wall velocity as well 
as the fusion effects of the nucleated bubbles. 

A Estimation of the Lepton Number Pro- 
duction Rate from the Bubble Wall 

Our starting Lagragian is that of Cohen et al. 0, namely 

£ = —C{standard model) + AC, (1) 

where 

AC = h^LtYd^NL - ^(PMxINl - H^lXdNI + {h.c). (2) 

Here V'l and Nl represent left-handed lepton doublets and right-handed 
neutrinos of G generations, respectively, stands for the additional singlet 
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Higgs-scalar, and H is the usual doublet one. The Majorana and Dirac mass 
terms of the neutrinos are given by the G x G mixing- matrices A]\^ and Xd, 
respectively. (C denotes the charge conjugation as usual.) 

Since the electroweak phase-transition is of first order, bubbles of the 
broken phase are nucleated within the unbroken vacuum when the phase 
transition begins and the system becomes supercooling. These nucleated 
bubbles grow and fuse with other bubbles; finally the whole space is filled 
up with the broken phase. Therefore we have the thermal non-equlibrium, 
one of the necessary conditions of the baryogenesis, during the development 
of this first order phase transition. 

We have, however, two kinds of Higgs scalar s, H and 0. 

If they acquire vacuum expectation values at different phases, then we 
have a complex configuration of the admixture of various kinds of bubbles. 
In order to avoid the complexity, we have assumed that there is only one 
kind of bubbles. For this purpose, the position-dependency of the vacuum 
expectation values {4>{x)) and near the interface of bubbles is asu- 

umed to be identical {4>{x)) oc {H^{x)). This rather unrealistic assumption 
is not too bad, if the lepton number production rate J'l comes mainly from 
the lighter components of neutrinos and not from the heavier ones the latter 
of which feel the configuration {4>{x)) sensitively. 

The Dirac equation of neutrinos for a given configuration of {4>{x)) and 
{H^{x)) reads from Eq(0) 

i{do - (T ■ V)^(x) - ia2M{x)^{xy = 0, (3) 

where 

vi>(x) = [u,,U2,---, i^G, N,, iV2, ■ ■ ■ , Ncf, (4) 

and (Tj(i = 1, 2, 3) are the usual Pauli matrices. 
In Eq.(^ the mass matrix M{x) is given by 

M(x)=(./' ^ , (5) 
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where ip{x) stands for the common configuration of {(f){x)) and {H^{x)) near 
the bubble wall, taking outside the bubble ( unbroken phase), 1 inside 
the bubble ( broken phase) , and the value in between within the wall of the 
bubble (interface of the two phases). Under the normalization adopted here, 

li^ohl < \{^M)kl\ (6) 

is understood for the see-saw mechanism to work. 

For a sufficiently large bubble of the radius R 3> with the typical 
temperature T ~ O{100GeV), we can reduce the three spatial variables to 
one, say z, representing the radial direction of the bubble. The decomposition 
oi^{z) in terms of the spin Sz, 



AEt 



r^{z)e^ 
tlji{z)e-'^^ 



+ 



^l)l{z)e 



-iEt 
iEt 



(7) 



indicates that the particle ■0i with — —1/2 is coming from the left, a part 
of which is refiected to the left as the anti-particle ^4 with the same spin of 
Sz = —1/2; Similarly the anti-particle 1^2 with Sz = +1/2 is coming from 
the left, a part of which is refiected as the particle ip^ with Sz = +1/2. 

We consider that the regions z < and z > 6^ belong to the unbroken 
and broken phases, respectively, while the region o < z < 5^, is occupied by 
the bubble wall. 

The refiection coefficients R and R for the particle {v, N) and its anti- 
particle {u, N) are defined by 

^^4(0) = RM^) and ^s{0) = ^^'2(0) (8) 

under the boundary comditions of V^i(oo) — and ^"2(00) = 0, respectively. 

In the derivation of R and R , we need to estimate the path ordered 
integral. 

I = Pexp (^i j^J dzA(ip(z) )^ (9) 
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with 2G X 2G constant matrix A. 
We approximate it to 



I A = exp (^i j^J dz^{z)A^ = e*''^, (10) 
where 77 is a parameter representing the shape of the wall {0 < rj < 1) ; 

1 fSw 

dzif{z) (11) 



Sw Jo 

It should be noted that the approximation (^) is exact when (f{z) changes 
linearly whitin the wall, < z < 6^, and for rj = 1/2. Therefore, the 
approximation (called linear one) is a rather good approximation. 

Now the reflection coefficiets of particle [{u, N) {u, N)] and anti-particle 
[{u, N) (z/, N)] read 

R ^ -U^DU and R = U^DU\ (12) 

where the diagonal 2G x 2G matrix D in the generation space denotes 

^ ^ ^ s + ^{sE - cr^Mf -^-')) x p-^ tan(^^P,, 

c - 2(cE - sr/Mf "^''"'^')) X P-i tan(5^P,) ■ 

Here the diagonal matrices c^, and are defined as 

= cosh 6'^, Sjj = sinh6'^, (14) 

with 

tanh 29^ = ^(Mf "^""'^'^ /E), (15) 

and 

Pr, = _ ^{diagonal) ^^^^ 
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where c = ci and s = si. The detail of the derivation of ( pTSD will be written 
elsewhere 0. 

Viewing the expression Eq.(|I^), we can understand that D is the reflection 
coefficient for each mass eigen-states {v^N) or {u,N), and it is sandwiched 
by the matrices ?7's which transform the weak eigen-states to the mass eigen- 
states, namely , 

UM{z in the broken phase)U^ = M^^°'^°"'°'^ . (17) 

Cooperation of the phase shift including in D and of the complex phases in 
the mass matrix (or in U) makes the difference between particlex (z/, A^) and 
anti-particles(i>, N). The number of the comples phases in U is G{G — 2) + 1, 
so that the minimum number of the generation is 2 for our purpose of the 
baryogenesis. 

We can read easily the following properties from our results Eg. ([13|) : 



i. |D(i?)p j{m^/Ey for m^/E 0, which means that for a massless 
neutrino the barrier or the wall disappears so that it perfectly transmit 
the barrier and does not be reflected. 

ii. \D[E)\'^ — > 1 for E < , that is, such neutrino is completely reflected 
by the "potential" barrier. 

iii. Near the threshold E = my + AE{AE > 0) , we have 



|B(£)|^.l_2,HM / + " (18) 



giving the cusp behavior \/E — m^, and damping rapidly for E > m^. 
This is the origin of "resonance behavior" in the baryon number pro- 
duction rate so called by Nelson et al. [Q. 
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iv. Asymptotically for E ^ rriy, 

\D{E)\' ~ {^)\\ cos' {6^E) + (i - r^)'sin2 {6^E)}, (19) 

where a series of "resonances" appear for S^E ^ 0. 

In general the phase transition occures at T = O{100GeV), which is 
approximately the physical Higgs mass mn, and wall width 6^ is roughly 
m]j^. Therefore, S^E, and are 0(1) and are not negligible, so that we 
will be troubled with the summing up a series of "resonances". 

Difference of the reflection rates between the two processes z/j — uj and 
z/i — i> Uj triggers the lepton number production 

A I p |2 I p |2 

= $:Jm(D,A*)x Jj^ (20) 

k^l 

where 

Jj^ ^ /m(^,,f/,,f/*f/*) (21) 

is the so called Jarlskog parameter [Q, typically representing the magnitude 
of CV violation arising from the mass matrix. The expression eq.(pU[) shows 
that in order for the CV violation to manifest itself, the other dynamical 
phase, Im{DkDi) is required to take part in the problem. 

Here we will meet with a difficult problem: The initial states i and the 
final states j are thermally averaged, so that if the thermal distribution is 
common for all i or for all j, then we have no lepton number production, 
= J2j ^ji = 0. This refiects the CVT invariance and the GIM can- 
cellation. The way to overcome this difficulty has been proposed by Farrar 
and Shaposhnikov 0. They has considered that the initial and final parti- 
cles have masses due to the finite temperature effects that lifts the common 
distribution for initial and final particles. 
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Here we restrict ourselves only to the process in which the neutrinos are 
coming from the unbroken phase and are reflected by the bubble wall. Then, 
initial and final neutrinos are massless at T = 0, but aquire the following 
finite temperature masses for and Ni, 

M.AT') = ^Az,AL (22) 

MNdT') = ^{XUmV, (23) 

(24) 

without breaking the chiral invariance [0]. 

Then, we have obtained the lepton number production rate fi per unit 
time generated from the unit area of the bubble wall, under the situation 
that the wall is moving toward the incoming neutrinos with the velocity v^j '■ 

1 r d^k^ 

h = - 7T;;TiE(l-/n(K))AL_(|A;f|)x^/^(E^), (25) 

with 

A-f'nm = ^(-^i ^ Li){\a\jRjiaim\^ — \a\jRjiaim\'^} (26) 

Here, Lj(j) denotes the lepton number of the neutrino un^j) ( or A^j^) ) in the 
flavor eigenstates, m{n) represent the mass eigenstates diagonalized by the 
matrix {aim}-, f{E) is the Fermi-Dirac distribution in the thermal frame, the 
energy and momentum with the affix W stand for the variables in the wall 
rest frame, and = (1 — ff,/c^)~^/^. 

Approximately we have the following expression 

. ^,^ Ml{T)-Ml{T) r d{kY/T) 

h ^ T — -j^{k^/T)ALU\k^\) 

X {In 7. + IH^JikYy + Ml{T)/T) +^E + \n2 

- 2^^^{kYy + Ml{T)/T} (27) 
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where 7^; is the Euler's constant (= —0.577- ■ ■). 

Then, Jl is the sum of 3 terms, having different dependence on the wall 
velocity (or 7^ ); 

h/T^ ^ (A In 7^ + - C7.) X J, (28) 

where the constants A, B and C depend on the model and the phase transi- 
tion temperature T . 

As an example, we take the 2 generation model, where only two lighter 
neutrinos Vi and V2 are assumed to contribute to the lepton number gener- 
ation. Using the parameters of mi = Mi(T) = lOOGeV", m2 = M2(T) = 
bOGeV and T = lOOGeV, we have 

A = 0.625 X 10"^ 
B = 0.110 xlO^^ 
C = 0.185 X lO-l 

(29) 

The same result is obtained in case of doubling the energy scale as mi = 
Mi(T) = 200GeV, m^ = M2(T) = lOOGeV and T = 200GeV. Diffficulty 
existing in the calculation of (|27|) is the integration over an infinite tower 
of small "resonances", so that the obtained values are preliminary. Using 
the thin wall approximation, there is no such trouble. Final lepton number 
production rate depends on J, the magnetude of the CV violation. 

B Temporal Development of the Phase Tran- 
sition Dynamics 

Next theme which we are going to study in this paper is the phase transition 
dynamics. During the development of the EW phase transiton, the bubbles of 
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the broken phase are nucleated at a rate of I per unit time and unit volume. 
They grow and fuse with each other. Finally the whole space is covered by 
the broken phase. The lepton number production occurs at the interfaces 
between the broken and unbroken phases, or the bubble walls. Therefore we 
need to know the temporal change of the total area A of these bubble walls, 
that is A = A{t). 

If the wall velocity is constant at any situation, then the total number 
of the lepton number production can be obtained without the knowledge of 

Ait); 

Nl= f v-^hM ■ v^A{t)dt = v-^fL{v^)Vtotai, (30) 



giving the lepton number density ri^ of the universe as 

n^L^=vZ'fL{v.). (31) 

The wall velocity v^^ is, however, by no means a constant, but is time- 
dependent 10; 

where R(t) is the radius of the bubble nucleated at time t, and Rc is the 
critical radius with which the bubble is nucleated. The F"^ is the friction 
coefficient which the global bubble feels when it grows in the heat bath; 

^^^^ 

where is the free energy of the Higgs field making the bubble. There- 
fore fa;(t) increases exponentially exp{2VT / R^) and approaches to the con- 
stant velocity of Vui{oc) = 2T/Rc. Then the typical length scale and time 
scale in this problem are Rc and to = i?^/2F. The F may be estimated using 
the linear response theory, and giving 0(T^^) ~ O(100Ge\^)~^). 
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The other important quantity is the nucleation rate / of the bubbles. We 
have the expression fl, 



I = T^{^^Y (34) 

with Fc{T) , the free energy of the bubble with the critical radius. Using the 
finite temperature effective potential estimated at 1-loop level reads 

V = ^(f)^ - ET(j)^ + D{T^ - T^)<p'^ (35) 

with 

D = -L(2m^ + ^| + 2m2) (36) 

E = ^^(2m3^ + 3m|) (37) 

T„ ~ (38) 

Aj. ~ A ~ -{mn/vY. (39) 

For the choice of values mw = 80GeV,mz = 90GeV,mt = 150(170)GeV^, 
V = 246GeV and the unknown parameter rriH = lOOGeV^ , we have D ~ 
0.27(0.33),E ~ 0.05, Tq ~ lOOGeV and A ~ 0.08. 

Correspondingly, we have roughly 

6^ ~ 2VX/ET ~ {9GeV)-^ (40) 
~ 2{DT)yX^/^-e (41) 

(42) 



and 
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with the latent heat 



(43) 



where is the temperature below which the broken phase appears with 
{ (j) ) = (pc- During the temperature Tc > T > Tq, the EW phase transition 
develops. The free energy of the critical bubble 

2^271 E^T^ 

gives the nucleation rate / as a function of x = Tc — T; Fig |I] depicts the func- 
tion / = I{x), where the vertical scale is normalized by R^to and x is in the 
unit of GeV. It is a problem of determining a fixed value of the temperature 
at which the first order phase transition develops; the problem should be an- 
swered by coupling the phase transition dynamics with the expansion of the 
universe. Here we try two typical value = lOOGeV and Tc2 = 200GeV for 
Tc. For these values, we have fixed the temperature at Ti = — 0.90GeV 
and T2 = Tc2 — 0.296*61^ , respectively, during the phase transition. 

We performed the computer simulation at Ti or T2 . In the simulations, 
we generated the critical bubbles at the rate of I(Ti or T2), these bubbles 
grow by changing their wall velocities according to (0), and the fusion effect 
of the bubbles is taken into account. 

The results of these simulation are given in Fig . and 0. The Fig ^ and 
1^ give the temporal evolution of the area of the wall A{t) and the volume 
V{t) of the broken phase, respectively. The value Ti and T2 are so chosen 
that the both simulations of A{t) and V{t) at Ti and T2 become identical. 

There is the exactly solvable model of Kolmogorov and Avrami 0, in 
which the v^^ and / are kept to be constant. This extremely attractive theory 
predicts for D = 3 as 

A{t)/V,otai = yt;S/ot^e-f^S/o*^ (45) 
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Our simulation differs considerably from the simple theory of the phase 
transition dynamics (|45|) . In our simulation, each portion i of the bubble 
walls gives different time-development, so that the finally produced lepton 
number density ul of our simulation becomes 

= E/ fL{vt{t)A\t)dt/Vtotai. (46) 

i 

Then, we have 

77^3 _ / —0.299 X 10~^ ■ J for the case 1 / 
~ I -0.303 X 10-2 . J for the case 2 , ^ ^ 

where the parameters used in the case 1 are Td = lOOGeV, Ti = Td — 
0.90GeV, mi = Mi(T) = lOOGeV, ma = M2(T) = 50GeV, and w^(oo) = 
0.482, whereas those used in the case 2 are = 200GeV, T2 = Tc2—0.29GeV 
, mi = Mi(T) = 200GeV , ma = M2(T) = lOOGeV , and v^{oo) = 0.311. 
The corresponding lepton number densities n° obtained in the Kolmogorov- 
Avrami model of the phase transition, read 

/m3 _ f —0.108 X 10^2 ■ J for the case 1 /.r.\ 

LI In onn i n-2 t f„„ 4-T,„ „„„„ o ^ ' 



-0.209 X 10^2 . J for the case 2 . 

Therefore, we have 

riL ( 2.769 for the case 1 



(491 

1.450 for the case 2 , ^ 



which is the most important results of our paper; The phase transition dy- 
namics with and without including the change of the wall velocity are found 
to affect considerably the final value of the lepton numer ( as well as the 
baryon number ) production. 

The detailed analysis will be given elsewhere O. 
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C Baryogenesis from Leptogenesis 



One of the difficulty in considering the baryogenesis at EW scale is the con- 
straint of mn < 4:5GeV which guarantees that the sphareron transiton is 
suppressed by the expansion of the universe and that the produced baryon 
number may not be washed away by the spharelon. 

In our case the sphaleron transition should be rapid enough to make the 
chemical equlibrium between the lepton number L and the baryon number 
B. Thanks to the B — L conservation, we have the non- vanishing equlibrium 
value for B. 

{B) = + 

= — -{1 + x) {produced L), (50) 

where (/a) = (Y) = is assumed with x = 0(1); x = —18/15 for the light 
Higgs doublet H and the light singlet scaler with three fermions of quarks 
and leptons and, x = —1/2 without this 0. Anyway we have ub at the 
same order of magnetude of the produced ul during the EW phase transition, 
without the unrealistic constraint of mn < 4:5GeV. This is the good point 
of this baryogenesis model induced by the neutrinos. To obtain a realistic 
value for ni/n.y , we need to have CP violation J of 0(10^^ ~ 10^'^). 
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Figure 1: nucleation rate I as a function of x — Tc — T [GeV]. The unit of 
the vertical axis isV — T^. 
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Figure 2: Temporal development of the total area of the wall for our simula- 
tion and for the Kolmgorov-Avrami theory (solid line: our simulation, dotted 
line: Kolmgorov-Avrami) 
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Figure 3: Temporal development of the volume fraction of the broken phase 
for our simulation and for the Kolmgorov-Avrami theory (solid line: simula- 
tion, dotted line: Kolmogorv-Avrami) 
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